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We apply new statistical tests on the four year 53 GHz 
DMR data and show that there is significant probability for 
the existence of coherent temperature discontinuities hidden 
in the CMB maps around the Galactic Poles. Comparing 
Monte Carlo simulations with the DMR maps we find that the 
probability for the existence of a coherent discontinuity with 



amplitude 0.2 x (- 



superposed on the gaussian fluctua- 



tions is more than double than the corresponding probability 
when no discontinuity is present. This result is consistent 
with the existence of a horizon size long cosmic string with 
mass per unit length /i given by G[iv s -y s ~ 10 -6 /7r where v s is 
the string velocity and 7 S is the corresponding Lorentz factor. 



I. INTRODUCTION 



New statistical tests have recently been proposed for 
the detection of coherent discontinuities (edges) hidden 
in CMB maps [Q . The main advantage of these statistics 
is that they focus on the large scale coherence properties 
of CMB maps and are therefore effective even in cases 
of low resolution maps provided that the sky area cov- 
ered is large. These properties are shared by the maps 
produced by the DMR instrument of COBE. The goal of 
this paper is the application of the new statistical tests 
analyzed in Ref. jl| on the four year DMR maps in the 
regions of the Galactic Poles (approximately 90° x 90° or 
32 x 32 pixel maps). The statistics calculated for these 
maps are the skewness, the kurtosis, the Sample Mean 
Difference (SMD) Q and the Maximum Sample Differ- 
ence (MSD) jl). These results are then compared with 
the corresponding results obtained from a large number 
of gaussianized DMR maps (similar power spectrum and 
random phases in Fourier space) with and without a co- 
herent discontinuity superposed. Thus for each applied 
statistic we find the value of the superposed coherent 
discontinuity amplitude 4jf = a x (4j£) rms that is most 
consistent with the actual DMR data. 

In particular we ask the following question: 'Assume 
that the DMR map around the North Galactic Pole is 
gaussian but has small coherent temperature discontinu- 
ities superposed on it. What is the most probable value 
of the amplitude of the superposed discontinuity? Also, 
what is the ratio of the probability that a discontinu- 
ity of the given amplitude is present over the probability 



that no discontinuity is present?'. Clearly, for a gaussian 
map we would find that the most probable value of the 
discontinuity amplitude is 0. 

Before describing the technique followed to address 
the above questions we will briefly review the statistical 
tests we have used (the notation used here is the same as 
in Ref. [Q). These tests involve both conventional statis- 
tics (skewness and kurtosis) and new statistics (MSD and 
SMD) optimized for the detection of coherent disconti- 
nuities in Id and 2d pixel maps. Assuming a 32 x 32 
standardized temperature pixel map TV, (i,j = 1, 32), 
the skewness s and the kurtosis k are defined as 



,s = <T 3 > = ]T^/(32) 2 

i,3 

fc = <T 4 >^]TT4/(32) 2 



(1) 
(2) 



h3 



These are conventional statistics and their values for 
gaussian maps with uncorrelated pixels are s = 0, k — 3. 
To define the statistics MSD and SMD we consider a par- 
tition of the CMB in two parts separated by a random 
line. In this study we have considered straight lines but 
the analysis can be generalized to more general types of 
lines with no significant changes in the results. Let k de- 
note the set of parameters that define the partition line 
and let T u and TJ be the mean temperatures of the two 
parts of the map (the indices l u' and T stand for 'upper' 
and 'lower' parts). The statistical variable Y±. is defined 



Tu-Ti 



(3) 



The statistics Sample Mean Difference (SMD) and Max- 
imum Sample Difference (MSD) are defined as : 



k 

MSD = max(y^) 



(4) 
(5) 



where N is the total number of partitions. Both MSD 
and SMD approach asymptotic values for large N. 

Consider now a gaussian pattern of temperature fluc- 
tuations with a small coherent temperature discontinuity 
defined by a partition k superposed on the map with co- 
herence scale comparable to the size of the pattern. In 
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Ref. pj| it was shown that the presence of the coherent 
discontinuity kg can be detected much more efficiently 
by the statistics SMD and MSD than by the skewness 
and kurtosis. A physically motivated mechanism which 
can lead to the production of a coherent discontinuity on 
CMB maps is the presence of a moving long cosmic string 
in our horizon [0. 

The main mechanism by which strings can produce 
CMB fluctuations on angular scales larger than 1° — 2° 
has been well studied both analytically || and simulated 
numerically [Q and is known as the Kaiser-Stebbins effect 
||. According to this effect, moving long strings present 
between the time of recombination t rec and the present 
time to, produce step-like temperature discontinuities be- 
tween photons that reach the observer through opposite 
sides of the string. These discontinuities are due to the 
peculiar nature of the spacetime around a long string 
which even though is locally flat, globally has the geom- 
etry of a cone with deficit angle 8irG^ (G is Newton's 
constant, \x is the mass per unit length of the string and 
we have used units with c = 1). The magnitude of the 
discontinuity is proportional to the deficit angle, to the 
string velocity v s and depends on the relative orientation 
between the unit vector along the string s and the unit 
photon wave- vector k. It is given by || : 

ST 

— = ±4irG[iv s j s k ■ (v s x s) (6) 

where "f s is the relativistic Lorentz factor and the sign 
changes when the string is crossed. The angular scale 
over which this discontinuity persists is given by the ra- 
dius of curvature of the string which according to simula- 
tions H is approximately equal to the horizon scale. The 
growth of the horizon from t rec to to results in a super- 
position of a large number of step-like temperature seeds 
of all sizes starting from about 2° (the angular size of the 
horizon at t rec ) to about 180° (the present horizon scale). 
By the Central Limit Theorem (CLT) this large number 
of superposed seeds results in a pattern of fluctuations 
that obeys gaussian statistics. Thus the probability dis- 
tribution for the temperature of each pixel of a CMB map 
with resolution larger than about 1° — 2° is a gaussian 

0- 

However, on large angular scales, despite the large 
number of superposed seeds there is also coherence of in- 
duced fluctuations. The fluctuations on these scales may 
be viewed as a superposition of a gaussian scale invariant 
background coming mainly from small scale seeds plus a 
small number of step-like discontinuities which are co- 
herent and persist on angular scales larger than 100°. 
The presence of coherent discontinuities due to late long 
strings, superposed on a gaussian background appears to 
be a generic feature of cosmic string models. In what fol- 
lows we will explore the consequences on non-gaussianity 
of this generic feature of the cosmic string models. 



II. STATISTICAL TESTS 

Most of the statistical tests for the detection of non- 
gaussianity focus on small scale properties (e.g. peaks) of 
CMB maps, because it is usually thought that the CLT 
tends to hide non-gaussianity on large scales due to the 
large number of superposed seeds on these scales. This 
expectation, however, is incorrect for the type of non- 
gaussianity induced by moving late long strings. The 
basic feature of this type of non-gaussian fluctuations is 
the large scale coherence and can thus be best detected 
by a specially designed statistical test. Such are the tests 
SMD and MSD analysed both analytically and numeri- 
cally in Ref. §. 

Here, we apply these statistics along with the more 
conventional skewness and kurtosis on the four year 53 
GHz DMR data. To find the probability that the DMR 
maps (around the Galactic Poles) are gaussian with a 
superposed discontinuity of a given amplitude 4jr = a x 
(4^) rms , we apply the following technique: 

1. We consider two 32 x 32 pixel maps of the 53 
GHz (A+B)/2 DMR data (monopole and dipole re- 
moved, optimum Galactic cut: 3881 surviving pix- 
els U) around the Galactic PolesQ (about 90° x 90° 
on the sky). We calculate the statistics: skewness, 
kurtosis, SMD and MSD. 

2. We construct a large number (1200) of gaussianized 
maps obtained from the DMR map by randomizing 
the phases in Fourier space and using a gaussian 
spectrum with mean and variance equal to the 
DMR measured spectrum. 

3. For each one of the constructed maps we calculate 
the four statistics: skewness, kurtosis, SMD and 
MSD. 

4. We superpose on the constucted maps a random 
coherent discontinuity of amplitude a x (4^) rms , 
obtained by a straight line partition and recalculate 
the skewness, kurtosis, SMD, and MSD statistics. 

5. Using the results from steps 3 and 4 we construct 
the probability distribution for the skewness, kur- 
tosis, SMD and MSD for various values of a (e.g. 
Figs. 1, 2). 



*As North and South Galactic Pole regions we have used 
faces and 5 of the COBE quadrilateralized spherical cube 
projection (CSC) in galactic pixelization, according to the 
plots and pixel numbering schemes found at: 

http:/ /www. gsfc.nasa.gov/astro/cobe/skymap_info. html. 
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6. We use the results from steps 1 and 5 to find 
the probability for obtaining the already calcu- 
lated, from step 1, COBE DMR values of skew- 
ness, kurtosis, SMD and MSD in the gaussian- 
ized maps as a function of a. Thus we have 
four probability distributions, one for each statistic: 
P s (a),P k (a),P S MD(a),PMSD(a) (Figs. 3, 4). 

For convenience we will hereafter use normalized 
probability distributions defined as: 

««> - « < 7 > 

Q(a) helps us answer the question: 'How probable is the 
existence of a coherent discontinuity within the studied 
maps?'. Fig. 1 shows the probability distribution of the 
skewness, for various values of a. Clearly, this statistic 
is insensitive to the detection of a stcplikc coherent dis- 
continuity in the data. The situation is similar for the 
kurtosis. 
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FIG. 1. Skewness probability distribution of the 1200 
gaussianized maps, obtained from the COBE DMR data 
around the North Galactic Pole, for two values of the dis- 
continuity amplitude a, (q=0 continuous line, a— 0.2 dashed 
line). The two dots represent the probability of getting the 
actual value of the skewness rendered by the DMR data for 
each value of a. 

Thus, these two conventional statistics cannot differ- 
entiate between a gaussian map and a map with a small 
superposed coherent discontinuity with a < 0.8. This is 
clearly seen in Figs. 3, 4 where Q s (a), Qfe(a) are shown 
to be quite flat functions of a. 

The situation is quite different for the MSD and SMD 
statistics. The probability distribution for the MSD, for 
example, is clearly very sensitive to the value of a and 
its maximum occurs for a = 0.2, thus favouring the ex- 
istence of a discontinuity (Fig. 2). 
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FIG. 2. MSD probability distribution around the North 
Galactic Pole (see caption in Fig.l). 

Clearly, the distributions Qmsd{o) and Qsmd{&), 
as it is shown in Figs. 3, 4, demonstrate a conspicuous 
peak at a ~ 0.2 with — 2.4 indicating that the 

presence of a discontinuity with a ~ 0.2 is more than 
twice as probable compared to a purely gaussian map. 
We therefore conclude that there are clear hints for the 
existence of a temperature discontinuity in the 53 GHz 
(A+B)/2 DMR data in the region of the North Galactic 
Pole superposed on the gaussian fluctuations (Fig. 3). 
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FIG. 3. Smoothed normalized probability distributions of 
the COBE DMR values calculated by the four statistics as a 
function of the discontinuity amplitude a. The distributions 
where obtained from the 1200 gaussianized maps covering the 
North Galactic Pole region. 

We have repeated the analysis for the data in the 
region of the South Galactic Pole (Fig. 4) with very 
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similar results, indicating the presence of a similar type 
of discontinuity. 
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FIG. 4. Smoothed normalized probability distributions, 
South Galactic Pole region (see caption in Fig. 3). 



discontinuities at the lcr level, but the probability for 
the existence of a coherent discontinuity with amplitude 
4? = a x (^jr)rms — 4 x 10~ 6 is more than double com- 
pared to the probability that no discontinuity is present. 
This result could be interpreted as an indication for the 
existence of a moving long string in our present horizon 
with GfiVs'js ~ lO -6 / 71 "- 
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The presence of such discontinuities on top of the an- 
ticipated gaussian fluctuations could be due to the pres- 
ence of a long string in our present horizon. According 
to the results from the application of the MSD statistic, 
such a string would go approximately through the galac- 
tic coordinates (230°. 98, 71°. 68) and (330°. 25, 48°.98) 
in the North Galactic Pole region (or pixels 612 and 
29 according to the quadrilateralized spherical cube). 
The corresponding discontinuity in the South Galactic 
Pole goes through the coordinates (315°. 00, —73°. 57) and 
(122°. 02, -73°. 12) or pixels 5733 and 5553 . It is straight- 
forward to check that the position of this discontinuity 
relative to the dipolc anisotropy excludes the possibility 
that the discontinuity is an artifact of improper removal 
of the dipole anisotropy. 



III. CONCLUSIONS 



We have shown that the 53 GHz (A+B)/2 DMR data 
in the regions around the Galactic Poles show clear hints 
for the existence of hidden temperature discontinuities 
coherent on large angular scales. These hints have been 
revealed by the application of the new MSD and SMD 
statistics which were especially designed for the detection 
of coherent discontinuities hidden in the fluctuation pat- 
terns. We have also shown that the conventional statis- 
tics of skewness and kurtosis are unable to reveal these 
features. The 53 GHz DMR maps in the regions around 
the Galactic Poles are consistent with the presence of no 
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